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Abstract 

We analyze the inverse problem, originally formulated by Dix [7] in geo- 
physics, of reconstructing the wave speed inside a domain from boundary 
measurements associated with the single scattering of seismic waves. We 
consider a domain M with a varying and possibly anisotropic wave speed 
which we model as a Riemannian metric g. For our data, we assume that M 
contains a dense set of point scatterers and that in a subset U C M, modeling 
a region containing measurement devices, we can measure the wave fronts of 
the single scattered waves diffracted from the point scatterers. The inverse 
problem we study is to recover the metric g in local coordinates anywhere 
on a set M C M up to an isometry (i.e. we recover the isometry type of 
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M). To do this we show that the shape operators related to wave fronts pro- 
duced by the point scatterers within M satisfy a certain system of differential 
equations which may be solved along geodesies of the metric. In this way, 
assuming we know g as well as the shape operator of the wave fronts in the 
region U, we may recover g in certain coordinate systems (e.g. Riemannian 
normal coordinates centered at point scatterers). This generalizes the well- 
known geophysical method of Dix to metrics which may depend on all spatial 
variables and be anisotropic. In particular, the novelty of this solution lies 
in the fact that it can be used to reconstruct the metric also in the presence 
of the caustics. 

Resume 

Nous analysons le probleme inverse, a l'origine formule par Dix [7] en 
geophysique, consistant a reconstruire la vitesse d'onde dans un domaine a 
partir des mesures aux frontieres associees kla. dispersion simple des ondes 
sismiques. Nous considerons un domaine M avec une vitesse d'onde vari- 
able et event uellement anisotrope que nous modelisons par une metrique 
Riemannienne g. Nous supposons que M contient une densite elevee de 
points diffractants et que dans un sous-ensemble U C M, correspondant a 
un domaine contenant les instruments de mesure, nous pouvons mesurer les 
fronts d'onde de la diffusion simple des ondes diffractees depuis les points 
diffractant. Le probleme inverse que nous etudions consiste a ^reconstru- 
ire la metrique g en coordonnees locales sur l'ensemble M C M modulo 
une isometrie (i.e. nous reconstruisons le type d'isometrie). Pour ce faire 
nous montrons que l'operateur de forme relatif aux fronts d'onde produits 
par les points diffractants dans M satisfait un certain systeme d'equations 
differentielles qui peut etre resolu le long des geodesiques de la metrique. De 
cette maniere, en supposant que nous connaissons g ainsi que l'operateur de 
forme des fronts d'onde dans la region U, nous pouvons retrouver g dans un 
certain systeme de coordonnees (e.g. coordonnees normales Riemannienne 
centrees aux points diffractant). Ceci generalise la methode geophysique de 
Dix a des metriques qui peuvent dependre de toutes les variables spatiales et 
etre anisotropes. En particulier, la nouveaute de cette solution est de pouvoir 
etre utilisee pour reconstruire la metrique, meme en presance des caustiques. 

Keywords: geometric inverse problems, Riemannian manifold, shape 
operator 
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1. Introduction: Motivation of the problem 

We consider a Riemannian manifold, (M,g), of dimension n with bound- 
ary DM. We analyze the inverse problem, originally formulated by Dix [7], 
aimed at reconstructing g from boundary measurements associated with the 
single scattering of seismic waves. When the waves produced by a source F 
are modeled by the solution of the wave equation (<9 t 2 — A g )u(x,t) = F(x,t) 
on (M, g), the geodesies 7^ on M correspond to the rays following the 
propagation of singularities by the parametrix corresponding with the wave 
operator on (M, g) and the metric distance d(xi, xq) of the points x\, 22 G M 
corresponds to the travel time of the waves from the point x\ to the point x%. 
The phase velocity in this case is given by v(x,a) = E" fc=1 g :jk (x)ajak] 1 ^ 2 , 
with a denoting the phase or cotangent direction. 

Below, we call the sets £ t)2/ = {j Vl v(t)', v G T y M, \\v\\ g = 1} generalized 
metric spheres (i.e. the images of the spheres {£ G T y M; ||£|| 9 = t) in the 
tangent space of radius t under the exponential map). The mathematical 
formulation of Dix's problem is then the following: Assume that 

1. We are given an open set T C dM, the metric tensor gjk\r of the 
boundary, and normal derivatives d?gjk\r for all p G Z + , where v is 
the normal vector of dM and g^ is the metric tensor in the boundary 
normal coordinates. 

2. For all x G T and t > 0, we are given at the point x the second 
fundamental form of the generalized metric sphere of (M, g) having the 
center y x>t and radius t. Here, y Xit = j X)V / x )(t) is the end point of the 
geodesic that starts from x in the ^-normal direction v(x) to dM and 
has the length t. 

Dix's inverse problem is the question if one can use these data to determine 
uniquely the metric tensor g on the set W C M that can be connected to V 
with a geodesic that does not intersect with a boundary. 

The above problem is the mathematical idealization of an imaging prob- 
lem encountered in geophysics where the goal is to determine the speed 
of waves (e.g. pressure waves) in a body from the external measurements. 
Roughly speaking, in Dix's inverse problem we assume that we observe wave 
fronts of the waves reflected from a large number of point scatterers, but for 
a given wave front we do not know from which scatterer it reflected. Let 
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us describe this problem next in more detail: Assume that the domain M 
that contains a quite dense set of point scatterers (called also diffraction 
points). Let T C dM be the acquisition surface on which we have sources 
and measure the scattered waves. Consider a point source at the point x G T 
that sends at time zero a wave that propagates into the domain and scatters 
from the point diffractors. In the single scattering approximation, all these 
point scatterers can be considered as new point sources. We observe on F 
the sum of the waves produced by these new point sources and the observed 
wave fronts coincide with the generalized metric spheres of (M, g). The wave 
fronts observed at the point x and the time 2t are produced by the point 
scatterers that can be connected to x with a geodesic whose length in the 
travel time metric is t. Combining the measurements corresponding to sev- 
eral point sources on T to simulate a measurement that would be obtained 
using a wave packet sent from T in the normal direction one can determine 
for given isT and t the second fundamental form (or equivalently, the shape 
operator) of the first wave front of the wave produced by the point scatterer 
at y Xt t that is observed at x to propagate to the direction — u(x). The wave 
speed in W should be then reconstructed using these shape operators. 

Also, we note that in some cases when in M there are surface discontinu- 
es instead of point scatterers, one can use signal processing of the measured 
data to produce the shape operators corresponding to point scatterers, see 
[13] . This recovery of point scatterer data, which we will not describe in 
detail here, involves differentiating with respect to the offset between source 
and receiver locations in dM. 

Earlier, Dix \7\ developed a procedure, with a formula, for reconstructing 
wave speed profiles in a half space 1R x M + with an isotropic metric that is 
one-dimensional (i.e. depends only on the depth coordinate). Despite this 
rather large restriction, Dix' algorithm has played a crucial role in imag- 
ing problems in Earth sciences. We generalize this approach to the case 
of multi-dimensional manifolds with general non-Euclidean metrics. Before 
continuing we add that since Dix, various adaptions have been considered 
to admit more general wavespeed functions in a half space. Some of these 
adaptions include the work of Shah [26], Hubral & Krey [IT] . Dubose, Jr. 
[8], Mann [20], and Iversen & Tygell [12]. In the case that direct travel times 
are measured, rather than quantities related to point diffractors, the related 
mathematical formulations are either the boundary or lens rigidity problems. 
Much work has been done on these problems (see [271 EE] and the references 
contained therein). 
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Figure 1: Notation used throughout the paper, following the geodesic 7. 

2. Introduction: Definitions and main results 

2.1. Background and notation 

Before continuing, we briefly mention some general references to Rieman- 
nian geometry (9j EH E5]- As this paper is intended also for researchers 
working on applied sciences, we recall some standard notations and con- 
structions in local coordinates, ( tic j '3C . . . . . OC ) . The metric tensor is given by 
gjk{x)dx^dx k and the inverse of the matrix [gj k ] is denoted by [g^ k \. Through- 
out the paper we use Einstein summation convention, summing over indexes 
that appear both as sub- and super-indexes. The Riemannian curvature 
tensor, Rijki, is given in coordinates by 

d d 

p« p« p« 1 rP p! pP -pi pP n pi p 

jkl ~ dx k jl dx l i k i l P k i k P h i kl ~ y ijkh 
where T l - k are the Christoffel symbols, 

3k 2 y \ dx k dxi Oxp J ' 

When X, Y e T X M are vectors, then the curvature operator R(X, Y) : 
T X M — y T X M is defined by the formula 

g(R(X, Y)V, W) = R ijkl X i Y j V k W l , V, W e T X M. 

Finally, = Va fc is the covariant derivative in the direction d k = gfs, which 
is defined for a (1,1) -tensor field Aj by 

d • 

and for a (l,0)-tensor field B l and a (0,l)-tensor field Bi by 

= ^- k B l + r[ p Bv, v kBl = A A _ V P k B p . 
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If / G C°° then the gradient of / with respect to g is a (l,0)-tensor field (i.e. 
a vector field) given in coordinates by 



(V/)' = g 



tj df 



2.2. Description of the problem and results 

Let us next formulate rigorously the setting for a modification of the 
above Dix's inverse problem that we will study in this paper. Let (M, g) 
be a (7^-smooth Riemannian manifold with boundary. We introduce an 
extension, (M,g), M C M, of (M,g) which is a complete or closed manifold 
containing M so that = g- For simplicity we simply write g = g and 
assume that we are given U = M \ M and the metric g on M \ M. We 
note that if M is compact and the boundary of M is convex, the travel times 
between boundary points determine the normal derivatives, of all orders, of 
the metric in boundary normal coordinates [15] ; hence, in the case of a convex 
boundary, the smooth extension can be constructed when the travel times 
between the boundary points are given. 

In the following we consider a complete or closed Riemannian manifold 
M with the measurement jlata givenjn an open subset U C M. The tangent 
and cotangent bundles of M at x G M are denoted by T X M and T*M, and the 
unit vectors at x are denoted by Q X M — {v G T X M; \\v\\ g = 1}. We denote 
by exp^. : T X M — > M the exponential map of (M, g) and when r] G fl x M, 
we denote the geodesies having the initial values (x, rf) by jx^it) = exp x (tv). 
Let Bj^(y, t) = {x G M; d(x, y) < t} denote the metric ball of radius t and 
center y and call the set {exp y (£) G M; ||£|| s = t} the generalized metric 
sphere. 

We now express the data described in the previous section in more pre- 
cise terms. Eventually it will be related to the shape operators of so-called 
spherical surfaces. 



Definition 1. Let U C M be an open set. The family Sff of oriented spher- 
ical surfaces is the set of all triples (t, S, v) satisfying the following properties 

(i) t > 0, 

(ii) S C U is a non-empty connected C°° -smooth (n — 1)- dimensional sub- 
manifold, 
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(Hi) There exists a y £ M and an open set Q C Q y M such that 

E = £j, )tj n = {l y , v (t)] rj £ fi}. (1) 

(iv) v is the unit normal vector field on E given by 

u(x) = jy lT ){t) at the point x = 7y lT) (i)- (2) 

The family Sjj of spherical surfaces is the set 

Su = {(t, E); there exists (t, E, v) £ S^}, 

that is, it contains the same information as Sff but not the orientation of 
the spherical surfaces. 

Note that if (t, E, v) £ S\J then in Q the set Q can be written in the form 
^ = {jx,v(x)(— t); x £ £} and as E is connected, thus also Q is connected. 

In reflection seismology one refers to E as the (partial) front of a point 
diffractor. Note that in the definition of Su we consider arbitrary y £ M, 
including points y in U . Due to this, we have the following result stating 
that Sjj determines both the metric g in U and the wave front data with 
orientation, that is, Sif. Even though the determination of the metric in U 
is not very interesting from the point of view of applications, we state the 
proposition for mathematical completeness. 

Proposition 2. Assume that we are given the open set U as a differentiable 
manifold and the family of spherical surfaces Srj. These data determine the 
metric g in U and the family of the oriented spherical surfaces S{J . 

Proposition [2] is provenjn the Appendix. 

Let x £ U and rj £ Q X M. We say that (t, E, v) £ S{J is associated to the 
pair (x, rj) if x £ E and v(x) = —rj. It is easy to see that if (i, E, u) £ S^f is 
associated to the pair (x, rj) then we can represent E in the form ([T]) where 
V = lx>q(t) and Q C Q y M is such that ( = —j X}V (t) £ Q. 

Once again, suppose that x £ U, rj £ Q X M. Now we proceed with more 
geometrical constructions along the geodesies j x>rj . Let F k {r) = F k (x,rj,r), 
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k = 1, 2, . . . , n be a linearly independent and parallel set of vector fields de- 
fined on 7 X)J? (R). This means that F k (x, 77, r) G T lx ^ r )M and v (r)F k (r) = 
0. We assume that F n (x,r),r) = %, v {r). Denote by gj k the inner products 

gjk = g{Fj(r),F k (r)). (3) 

Because the vector fields Fj are parallel, ^-fe does not depend on r. Let 
P = /-'(x, 77, r), j = 1, . . . , n be the co-frame dual to Fj. This means that 



(f{r),F k (r)) = & 



where (., .) denotes the usual pairing of T y M and T*M . Let \& = ^ x ,ij '■ R n — >■ 
M be the map 



'n-l 



^(s 1 ^ 2 ,...^ 71 \r) =exp ?(r) ^VF fc (x,77,r) , where g(r) = 7^0)- 



For all r G R the point g(r) has a neighborhood Bj^(q(r), e), e > so that 
there exists an smooth inverse map : B^(q(r),e) —> R n . We call such 
inverse maps the Fermi coordinates. 

Let i? be the curvature operator of (M,g). Below, we denote 

0, 7i, r) = (f k (x } rj, r), R(Fj(x, T], r), %,r,(r))% tri (r)) (4) 

and call r^(x,^,r) the curvature coefficients of the frame (Fj(x, rj, r))" =1 . 

For t > 0, let C(x,r],t) be the set of those r > for which 2 r = T^r) 
and 2/j = ■y x ,r 1 (t) are conjugate points on the geodesic 7 Xjr? . When t > r > 
and r ^ C(x, 77, t), there is a spherical surface £ r ^ := £ t _ r)J/ti n of the form (jl| 
with some neighborhood Q C Q yt M of £ = — T^i). Then z r = 7 x ,r;( r ) G 
E r>t . In this case we write the shape operator of S r t at z r as S XiVtr>t . Thus 
S x ,r),r,t G (Tl) Zr M is defined by 

S x ,r),r,tX = Vjf 

for all X G T Zr M where z/ is the normal vector field for S rt satisfying (2J). 
Let us write the shape operator with respect to the parallel frame: 

S x , v ,r,t = sJ(x,?7,r,t)/ J (x,r7,r) <g> F k (x,r),r). (5) 
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We say that Sj(x, r], r, t) are the coefficients of the second fundamental forms 
of the spherical surfaces on the geodesic corresponding to the frame 
Fk(0). The family of the oriented spherical surfaces Sy and the metric tensor 
g in U (which are in fact determined by Su by Proposition [2])), determine 
all triples (t, E, v) G S^f that are associated to the pair (x,rj). Note that 
if t > is such that x and r y XtV {t) are not conjugate along j x>T) , there exists 
at least one triple (t, E, v) G S^J that is associated to (x,rj). Note that all 
such surface E have the same the shape operator S XtV ^ t with r = at x. 
Thus, by computing the shape operator of such surfaces E at x we can find 
the operator S XiVt o it and furthermore the coefficients s^(x, r], 0, t) for all t > 
such that 7x,r?(^) is not conjugate point to x. 
Our main result is the following. 

Theorem 3. Let (M, g) be a complete or closed Riemannian manifold of 
dimension n and U C M be open. Then 

(i) Let x G U , r\ G Q X M and F k {r), k = 1, 2, . . . , n be linearly independent 
parallel vector fields along j x ^. Assume that we are given gjk = g(Fj, F k ) for 
j, k = 1, 2, . . . , n and the coefficients of the second fundamental forms of the 
spherical surfaces corresponding to the frame Fk(0), that is, Sj (x, 77, 0, t), for 
all t G K + \ C(x, 77,0). Then we can determine uniquely Sj(x,r],r,t) for all 
t > and r < t, r G M + \ C(x,r],t), and the curvature coefficients r^(x, r],r) 
for all r G M + . 

Consequently, 

(ii) Assume that we are given the open set U C M , metric g on U , Xq G U 
and a unit vector rjo G fl XQ M and let V be a neighborhood of (xo,r] Q ) in TM. 
Moreover, assume we are given the set 

S° r v := {(t, E, v) G S° r ; (x, -u(x)) G V for all x G E}. 

Then for all r > there is p = p(xQ,r]Q,r) such that the Fermi coordinates 
^xo,vo associated t° the geodesic 7 XOjW (IR) are well defined in an open set 

V X0 , V0 ,r = y xo , m {B R n-i(0,p) x (r- p,r + p)), 

and the above data determine uniquely (^> XOtVO )^g, that is, the metric g in 
Fermi coordinates in V XOtTt0tr . This is the meaning of "reconstruction of the 
isometry type of the metric in the Fermi coordinates near 7 XOir)0 (]R + ). " 
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In the above theorem, (i) says that the shape operators S x ^^ r>t of the 
lx,r)(t) centered generalized spheres with radius t — r can be uniquely de- 
termined from the shape operators S X)V .oj corresponding to r = when the 
metric in U is known. The claim (ii) says that the Riemannian metric near 
the geodesic ^ XOtm {K + ) can be determined from the knowledge of wave fronts 
propagating close to this geodesic. From the point of view of applications, it 
is particularly important that the reconstruction can be done past the first 
conjugate point, that is, beyond the caustics of the reflected waves. 

We point out that the reconstruction method we develop in this paper 
is constructive and that it is based on solving a system of ordinary differen- 
tial equations which are satisfied by dfSj(x,T),r,t), p — 0, 1,2,3 along each 
geodesic. 

Using Theorem [3] we can prove the unique determination of the universal 
covering space. 

Theorem 4. Let M and M' be two smooth (compact or complete) Rieman- 
nian manifolds and U C M and U' C M' be such non-empty open sets that 
there is a diffeomorphism $ : U — > U' and Su> = {(£, $(E)); (t, E) G Su}. 
Then there is a Riemannian manifold (N,gw) such that there are Rieman- 
nian covering maps F : N — >■ M and F' : N — )■ M' , that is, M and M' have 
isometric universal covering spaces. 

The next example shows that Su does not determine the manifold (M, g) 
but only its universal covering space. 

Example 1. Let (M, g) be the flat torus R 2 /Z 2 , and (M',g) be the 
flat torus R 2 /(2Z x 2Z). Below, we consider (M,g) and {M',g') as the 
squares [0,1] 2 and [0,2] 2 with the parallel sides being glued together. Both 
(M,g) and (M',g') have the universal covering space M 2 with the Euclidean 
metric. Let U and U' be the disc B(p, |) of the radius | and the center 

p = (g, |). We see that both the collection for (M,g) and the collection 

Sfj, for (M',g f ) consist of triples (E, t, v) where t > 0, E is a connected 
circular arc having radius t that is a subset of the disc B(p, |), that is, 
E = {(tsina + x, tcosa + y) G U\ \a — a \ < c }, and v is the exterior 
normal vector of E. Note that the circular arc E can also be a whole circle if 
t is small enough. This shows that the knowledge of U and S^, is not always 
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enough to determine uniquely the manifold (M, g) but only its universal 
covering space. 

The above example shows that S^J, contains less information than many 
other data sets used in inverse problems that determine the manifold (M, g) 
uniquely such as the hyperbolic Dirichlet-to-Neumann map A wave , the parabolic 
Dirichlet-to-Neumann map associated to the heat kernel considered in jH [16l 
IT71 fT8l l2T| |29| [30] . the boundary distance representation R(M) (i.e., the 
boundary distance functions) considered in [3l [191 EZ] > or the broken scatter- 
ing relation considered in [TJ]; see also related data sets in [32] . 

2.3. Jacobi and Riccati equations 

Before moving to the actual reconstruction procedure we collect a few 
more geometrical formulae that will be useful. First, if we fix the initial data 
(x, rj) for the geodesic and a t > 0, then S t (r) = S XiV>r; t can be thought of as 
a (l,l)-tensor field on the geodesic y X)V . Let y t = y x>v (t), z r = lx,rj( r ), and 
( = —(t — r)7 Xj?? (t) so that z r = exp yt ((). Assume that y t and z r are not 

conjugate points along y Xtn {t). Then the exponential function exp^ t : T yt M — > 

M has a local inverse F t = exp" 1 in a neighborhood V of z r and the function 
ft{z) = \\F t (z)\\ g is a generalized distance function (i.e. ||V/t(^)|| g = 1, 
z G V). As the spherical surface Yi y j- r ,n near z r can be written as a level 
set of a generalized distance function f t , it follows from the radial curvature 
(Riccati) equation [211 sect. 4.2, Thm 2] that 

-V dr S t (r) + S t (r) 2 = -R dr (r), (6) 

where R dr {r) : T Z M -> T Zr M, R dr (r) : V \-> R{V,d r )d r is the so-called 
directional curvature operator associated with the Riemannian curvature R 



of (M, g) and d r = —V ft- Note that at the point 7^ 



r) we have d r = Tz^r). 



Let Tj(r) = rj(x,r),r) be the coefficients defined in (ij), that is, 

v){r) = {f\R{F v d r )d r ). 

If we also express S t (r) = S x „ r j in the parallel frame as in ^ with s*(r, t) = 
Sj(x,rj,r,t), then equation (|6[) becomes 

-^sj + sjs; = -rj. (7) 



12 



Also there is an equation we will use relating Jacobi fields along 7^ to S t . 
A Jacobi field J(r) along the geodesic 7^ is a vector field satisfying 

V 2 dr J + R(J,d r )d r = 0. (8) 

Writing 

J = J k F k 

this equation is 

d 2 r J k + rp j = 0. 

Finally, it follows from [TUl p. 36] that if J\ r=t = 0, then V_a r J = S t J, that 
is, 

-W dr J = S t J. (9) 



With respect to the parallel frame equation (|9j) reads as 

- 3 r J k = Sp j . (10) 

Our strategy for reconstruction is to show that from the data we can 
reconstruct shape operator S X)rjtrj t along each j XjTI using the Riccati equation 
^ . Then the Jacobi fields may be calculated using (J8]) , and since the Jacobi 
fields are the coordinate vectors for the local coordinates introduced in the 
next section this then allows recovery of the metric with respect to those 
coordinates. 



3. Coordinates associated to a spherical surface S 

We now introduce a set of coordinates in which we will initially jrecover 
the metric g. Suppose that we fix x$ G U — M\M and 770 £ T Xo M, and 
then pick a large to > such that xo and 7 a;o , ??0 (£o) are n °t conjugate points 
along 7x , r? - 

We will use the notation yt — O^o^o (^o)- Suppose that So C U 
is a spherical surface containing xq given by y to and t in the form Q. When 
v = u(x) is the normal vector field of S oriented so that v(xq) = —i]q, we 
have (£ ,t ,v) e Stf 

Let us take arbitrary coordinates X : V — >■ IR n_1 on V C So such that V 
is a neighborhood of xq and X(xq) = 0. To slightly simplify the notation, 
we identify V with its image X(V) C M n_1 and the points x & V with their 
coordinates x = (x 1 , . . . ,x n ~ 1 ) = X(x). Also, denote the geodesic starting 
normally to S from X~ l (x) G E by Jx(t) = 1x-u(x)(t). Let 1 : S — > M be 
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the identical embedding. For x G X(V), we define at the point q = X 1 (x) 
the vectors 

Fj(x) = l* ( -4U J , j < n - 1, and F n (x) = -z/(X _1 (x)) 
V dXi J 

Then (Fj(x))] = l is a basis for T g M. Let Fj(x,r) be the parallel translation 
of the vector Fj(x) along the geodesic 7$(r). 

Let be the set of those r G M + for which x G E and 7$(r) are 

conjugate points on the geodesic 73, and let 

W = {(£, r) G V x M + ; ifF, t G M+ \ C(£)}, 
W = { 7$ (r) G M; (x,r) G W}. 

Then the map 

X 9 :W-*WCM, X 9 (x,r)= 1 z(r) (11) 

is a local diffeomorphism. Below, we use the local inverse maps of X^ as local 
coordinaes on W. The (x, r) coordinates, basically, are Riemannian normal 
coordinates centered at yt , but parametrized in a particular way: x can be 
thought of as a parametrization of part of the sphere of radius to i n T ytQ M, 

and then r corresponds to the radial variable in T yto M. Note also that the 
coordinate vectors in these coordinates are Jacobi fields along the geodesies 

lx- 

Below, we also denote 

sj (x, r,t)=s k J (X~\x),v(X~ 1 (x) ) , r, t) , 
v){x,r) =v k ] {X~\x)^{X-\x)) 1 r). 



4. Reconstruction of the shape operator along one geodesic 

Next, we fix x and aim to reconstruct the shape operator along 7$. 
Throughout the section we will suppress the dependence of all quantities 
on x because it is considered to be fixed. 

To simplify notations, in this section we use the conventions 7 = 7$ and 

Y k .(r-) = r k (x,r). 
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4-1. A lemma concerning curvature 

Our first step toward the reconstruction will be to prove a lemma relating, 
roughly speaking, the inverse of the shape operator of spherical surfaces as 
their radius goes to zero with the directional curvature operator at the center 
of the spherical surfaces. 

Lemma 5. Let S(r, t) = (s^,(r, t))™^ be given by the matrices defined in 

ti > and i\ be the injectivity radius of {M,g) at 7(^1). Let t,r G 
\t\ — ii/2, ti] with t > r and K(r, t) = S(r, i) -1 . Then 

K(r, t) = (t-r)I+ R(t) + 0((t - r) 4 ), R(t) = (rj(*) (12) 

where 0((t — r) 4 ) is estimated in a norm on the space of matrices R nxn . 



Proof. Throughout this proof all the indices and sums will run from 1 to 
n — 1 unless otherwise noted. We use Jacobi fields j k m \(s; r, t)F k (s) on 7([r, £]) 
satisfying 

^jf m) («;r,t)+rJ( a )jf ro) ( S ;r l t) = 0, s G [r,f], (13) 
supplemented with the boundary data 

j k [m) (s;r,t)\ s=r = 5^ i k {m) (s;r,t)\ s=t = 0. (14) 

We will consider these when i — r = e > is sufficiently small. Next, we 
freeze t and introduce notations 



>>*) = j{ m )(i- z ;i- £ ^). 



= >i(t-z)> 



where 2 G [0, e]. Equations (13)-(14) attain the form 

^«(m)(2.*)+^*K m) («.*) = 0, *e[0,e], 

v\ m) (z,t)\ z=0 = 0, uf m) (^,*)| a 

We then let 

w k m) . e {y,t) = v k m) (ey,t), 

vlAvit) = pL(ey,t), 



(15) 
(16) 

(17) 



(18) 
(19) 
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where y G [0,1]. We drop the subscripts and superscripts for simplicity of 
notation and view w and a as matrices. Then 

d 2 w(y,t)+e 2 a(y,t)w(y,t) = 0, ye [0,1], (20) 
w(y,t)\ y=0 = 0, w(y,t)\ y=1 = I. 

The supremum of the norm of the Riemannian curvature tensor is bounded 
in compact sets and thus we see that ||cr(y,i)|| is uniformly bounded over 
y G [0, 1] and t G [t — i /2,t ]. Thus we see that there are C\ and e 1 > 
such that if < £ < £i, then 



+ e 2 a(y,t)) i^ao,!])^ 1 ^,!]) < Ci, (21) 



for some constant C\ > 0. This implies that there is a C2 > such that for 
all t and e corresponding with r, t G [to — «o/2, to]) 

|K-,t)||tfi([o,i]) <C 2 . (22) 

We expand a(y,t), 

a(y,t) = a(0,t)+£ a (y,t), (23) 

with 

HM-,*)IU~([o,i]) <Cs£, (24) 
where C3 depends on the supremum of ||V-R|| ff and the norms of f k and Fk 
on the geodesic 7([0, t ]), i.e., the C 3 -norm of the metric. We expand w 
accordingly, 

w(y,t) = w°(y,t)+£ w (y,t), (25) 

where 

d 2 y w°(y,t)+e 2 a(0,t)w°(y,t) = 0, ye [0,1], (26) 

f)|»=o = 0, w°(y,t)\ y=1 = L 

We observe that there is a constant C 4 > 0, such that for all y G [0, 1] and 
t e [*i-ii/2,fi], 

||(^ + e 2 tr(0, *)) _1 Hi»([o,i])^ao,i]) ^ ^- (27) 
This implies that there is a C5 > for all t and e such that r, t G [ii — ii/2, ti], 

lk°(-,*)ll^(P,i])<^5. 



16 



Now 



d 2 £ w (y,t)+e 2 a(0,t)£ w (y,t) = -e 2 £ a (y,t)w(y,t), ye [0,1], 
£w(y,t)\v=o = 0, £ w (y,t)\ y=1 = 0. 



Using (24), (22) and (27), we find that there is a C 6 > such that 
for e sufficiently small. 



Denoting A = A(t) = a/o"(0, t) e C5 n 1 ) x ( n x ) (we can use any branch of 
the matrix square root) we get 

w°{y,t) = [sin(eA(t))]- 1 sin(eA(Oy). 
Expanding this solution in e yields 

w°(y, t) = y ( I + ie 2 A(t) a (l - y 2 ) + 0(^ 4 ) ) 

whence 

t%, = y ( / + ^ 2 A(t) 2 (l - y 2 ) ) + £ w;1 (y, t), (28) 
where , t) ||hi([o,i]) < C' 6 e 3 , and 



d y w(y) = 1+ ^e 2 \(t) 2 - l -e 2 \{t) 2 y 2 + d y £ W]1 (y, t), 



where |]9 v £ tu; i(- , t) ||£a([o,i]) < CVr 3 - 
We recall (15) and differentiate, 



d z v k m) (z,t)\ z=£ = -d s j k m) (s;t- e,t)\ s=t _ £ . 
Because v k (m) (z, t) \ z=0 = jf m) (t; t - e, t) = 0, 

-d s 3 k m) (s;t-e,t) = s k (s,t)f (m) (s;t- e,t) 



(cf. (10)). Using this identity at s = t — e, we find that 

d z v(z, t)\ z=£ = S(t-e, t)v(z, t)\ z=£ . 



(29) 



(30) 
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However, v(z,t)\ z=£ = I, so that, with S(t - e,*) := {gikS k m {t - M))z)mLi> 

S(t-e,f) = S(t-e,tMM)-«(M) = "KM) 
{d z v{z, t) ■ d z v(z, t) + d z v(z, t) ■ v(z, t)) dz, 



,k „J 



where v-v stands for v^v^dkj- Substituting (18) yields 
&(t-e,t) = I (e' 1 d y w(y,t)-e' 1 d y w(y,t)+e~ 2 dyw(y,t)-w(y,t))edy 



o 



s 1 (dyw{y,t)-d y w{y,t) - e 2 a{y,t)w{y,t) -w{y,t)) dy, 
Jo 



using Jacobi equation (20). 



Inserting expansions (28) and (23) gives 



1 



S(t-e,t) = s- 1 I ((/+^ 2 A(t) 2 -^ 2 A(t)VJ" (I+^\(ty-^\(t)V 

J 

- 5V(o,t)y 2 (/ + \e 2 \mi - y 2 )} (/ + \z 2 Kt)V - y 2 )))dy + 0(e 2 ) 



so that 



S(t-e,t) = e 



Then, as p(0,t) = a(0,t), we obtain 



so that 



S(t - s, t)- 1 = e ( I + j P {0, t) + 0(e 3 



K(r, t) = (t-r)I+ (t / )3 R(t) + 0((t - r) 4 ) 



using that p(Q,t) = R(i). 
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4-2. Reconstruction 

In this section we complete the proof of Theorem [3} The main part is the 
proof of claim (i) which is given in the following proposition. 

Proposition 6. Functions s k (0,t) , t > 0, determine uniquely functions r k (r) 
for r G [0, to] ; an d s k (r,t) for r,t e [0,t ] with r <t where it is defined. 



Proof. We are given the matrices S(0,t) = (s£(0, *))"fc=i> t > 0. Using 
Lemma 5, it follows that the curvature matrix, R(r) = (r^r))™^, satisfies 

R(r) = ^ t 3 K(r,t)| t=r ; (31) 

similarly, R(r) = — ^<9 T 3 K(r, t)\ r=t . 
Using 

d r S(r,t) = S(r,t) 2 + R(r) 

(cf. Q) we find that 

d r K(r,t) = -(S(r,t))- 1 9 r S(r,t)(S(r,t))- 1 

= -(S(r,t))- 1 (S(r,t) 2 + R(r))(S(r,t))- 1 

= -J-(S(r,t))- 1 R(r)(S(r,t))- 1 

= -J-K(r,t)R(r)K(r,f). 

We let dt act on the final equation above, and obtain 

d r ((d t K)(r,t)) = d t (-I-K(r,t)R(r)K(r,t)) 

= -((d t K)(r, t)R(r)K(r, t) + K(r, *)R(r)(d t K)(r, t)). 

Computing the second and third i-derivatives in a similar manner, and de- 
noting V = V(r,t) = <y j (r,t)) 3 =0 , V j (r,t) = d{~K(r,t) and R = R(r), we 
obtain the equations 

d r V° = -I-V°RV° } (32) 

d r V l = -(V^V + V°RV 1 ), (33) 

d r V 2 = -(V r2 RU° + U°RV r2 + 2U 1 RU 1 ), (34) 

d r V 3 = -(V 3 RV° + V°RV 3 + 3V 2 ~RV 1 + 3U 1 RU 2 ). (35) 
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Since R depends on V 3 , this system is "closed". We define the operator T 
by 

(7V)(r) = V 3 (r,r) (36) 



so that, with (J3TJ) , R = R(r) = |(TV)(r). Hence, 

d r V° = -I - ^V°(TV)V°, (37) 
drV 1 = -hv 1 (TV)V° + V\TV)V l ), 
d r V 2 = -^(V 2 (TV)V° + V°(TV)V 2 + 2V 1 (TV)V 1 ), 
d r V 3 = -^(V 3 (TV)V° + V\TV)V 3 + 3V 2 (TV)V 1 + 3V\TV)V 2 ). 
We write this as 

d r V(r,t) = F(V(r,t),(TV)(r)), 
where the map F is a polynomial of its variables. We then introduce 

7: W(r,t) h-> F(W(r,t),(TW)(r)) 



so that the system (37) of nonlinear differential equations attains the form 

d r V(r,t) = (JV)(r,t). (38) 
Assuming that we are given S(0, t) with t > 0, we know the initial data 

V (t) = V(0,t) = (diiS^t))- 1 ) 3 (39) 



We now address whether the initial value problem (38)-(39) has a unique 
solution. 

Let us now assume that t\ > is such that the geodesic 7([0,ti]) has 
no focal points. This implies that the matrix K(r, t) is bounded on (r, t) G 
[0, ti} 2 . Let /C be a prior bound on the Riemannian curvature, that is, < 
K, on 7([0, ti]). The constant K. depends thus on the C 3 -norm of the metric. 
Using Cor. 2.4], which implies that 



S(r,t)> CQs(v ^ (t - r)) ^J 
sin(^(t - r)) 
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boundedness is then guaranteed if t — r < 7r/(4v/C). Let < t 2 < t\ and 

Y t2 = c([oMr-,c([oM^ (n ' 1)x{n - 1) ) 4 ) 

equipped with the norm 

\\ V \W t2 ■= SUp \\V(r, •)||c7([0,t2];M(»- 1 )x("-i)) 4 
re[0,t 2 ] 

= sup max ||y 3 (r,t)|| K (n-i)x(n-i). 
(r-,t)e[o i t 2 ] 2 - 7 ' e 't '--' 3 ^ 

It is immediate that 



V 6 (r,r)\< sup ||yV,t)|| R („-i)x (n -i) < 

(r,t)e[0,i 2 ] 2 



'2 ' 



If B t2 (1Z) C Yj 2 is the zero centered ball of radius 7£ > 1 in Y t2 , because T 
contains no differentiation, we find that 



T: B t2 (1Z) -> Y t 



t 2 



is (locally) Lipschitz, with Lipschitz constant L(TZ), that does not depend on 
We reformulate the differential equations in integral form, HV = V, with 

rr 

H : Y t2 ^Y t2 , {HW){r,t) = V {t) + / F{W{r\ t))dr', r,te[0,t 2 ]. 

J o 

Clearly, H : B t2 (1Z) — > Y t2 is (locally) Lipschitz, with Lipschitz constant 
t 2 L(1Z). For if to be a contraction, we need that 

t 2 L(K) < 1. 

To guarantee that H(B t2 (7Z)) C B t2 (lZ), we require that 

ll^ / o||c([0,ii];R(™- 1 ) x ("- 1 )) 4 +*2 (1 + 4:1Z 3 ) < 1Z. 

We choose 

1Z = 2 || Vo||cQo,ti];R( n - 1 ) x (' i - 1 )) 4 + 1 

where the norm || Vo 1 1 c?([o,ti.] -rc™— 13 x c«— 1) )■* can al so be bounded in terms of /C 
using Gronwall's lemma. Indeed, using (2U Cor. 2.4], equations (32)-(35), 
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Lemma [5j and Gronwall' s lemma, we can obtain the following basic "for- 
ward" estimates: 

\\V (r,t)\\ R ( n -i)x (n -i) ^ 
II^^OIkCn-DxCn-D <e Tl/3 , 

\\V 2 (r,t)\\ R(n - 1)x(n - 1) <21Ce 4T \ 

and 

||V 3 (r,t)|| R( n-i)x(n-i) < 8/C 3 (l + ne 8ri )e Tl 

for r,i 6 [0, ii] and t — r < T\, t\ = Ti(/C) = 7r/(4vX); the maximum of 
these results is an estimate for 1 1 Vq 1 1 cjrro.tiljMC 1 *- 1 ) x C"- 1 ))'* m terms of K. Then 
we choose 

1 / 71 1 71 \ 

t2 = 2 mm ^'w^rr^)J' (40) 

we see using the Banach fixed point theorem that H has a unique fixed point 



in B t2 (JZ). Thus (38)-(39) has a unique solution V G Yt 2 . 

Let C r be the set of those t G [0, oo) for which j(r) and 7(t) are conjugate 
points. Recall that we assume that we are given the functions s* (0,t) for all 
t G R + \ C and Z, fc G {0, 1, 2, . . . , n}. Thus we know S(0, t) and K(0, t) = 



S(0,t) _1 , for all t6E + \C . Let us choose t 2 > so that (40) is valid. Then 



the equation (38) has a solution V(r, t), (r, t) G A t2 := {(r, t) G [0, t 2 ] , r < t} 



with initial data V(0, t) = (9jK(0, t)) 3 =0 . Solving for V(r,t) using equation 
(38) gives us also the curvature matrix R(r) for r G [0, t 2 ] by applying T 
(i.e. we compute R = \TV). We will now switch notation replacing t 2 by t\ 
intuitively indicating that for the first step we can reconstruct R from to 

h. 

Next we do step by step reconstruction showing how to reconstruct R on 
the whole interval from to to i n steps. First, we observe that for t G M. + 
outside of the discrete set C r for which 7(7") and j(t) are conjugate points, 
the matrix S(r, t) is well defined. 

For given t > t\, t ^ Co, we will next reconstruct S(r, t), r G [0,ti]. 
As we know R(r) for r G [0,£i] and the matrices S(0, t), we find on the 
interval r G [0, ti] the solutions j fc (r, t) of the Cauchy problems for the Jacobi 
equations 

9 r 2 jf(r,t) + rJ(r)jf(r,t) = 0, r e [0,h], (41) 
j?(r,t)\ r=0 = 6f, 9 r jf(r,*)|r=o = -sf(0,*). 
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Now, when t is given, for all r G [0, ti] \ C t the vectors {jf (r, t)}™ =1 are linearly 



independent; then, the equations — d r jj(r, t) = s 3 k (r, t)jf (r, t) (cf. (10)) deter- 
mine S(r, t). Summarizing the above, for all t G IR+ \ Co and r G [0, t\) \ C t 
we can determine S(r, t). As t S(r, t) is continuous for t G M + \ C r , we see 
that we can find S(r,t) for all r G [0,ti] and t > r such that t G M + \ C r . 
In particular we can determine S(ti,t) for all t > ti such that t G R + \ C tl . 
This yields a new dataset in the interior of M at the point 7(ti). We now 
repeat the above argument with replaced by t\ to recover R and S(r, t) on 
another interval [ti,^]- 

As the metric is assumed to be C 3 -smooth, the size of the steps (i.e. 
t2 — ti) are in a compact set uniformly bounded below by the right hand 



side of (40). Thus we can complete the reconstruction in a finite number of 
steps and the proof is complete. This completes the proof of the claim (i) of 
Theorem HI 

5. Reconstruction of the metric in Fermi coordinates and the re- 
construction the universal covering space 



We are now ready to prove claim (ii) of Theorem |3} 

Proof. Let us now recall jsome considerations from section [3] For any xq G 
U = M \ M and rj G T Xo M and t > we considered the spherical surface 
So C U with a center y to = 7x ,r) (^o)- On E we considered the coordinates 
X : V — > IR n_1 in a neighborhood V C So of xo and on the geodesies 7a(r) 
we defined the parallel frames Fj(x, r), j = 1, 2, . . . , n. Note that as S C U 
and we know the metric tensor g on U, we can determine the inner products 

g jk {x)=g{F J {x),F k {x)), (42) 

and g(Fj(x, r), F k (x, r)) = gj k (x) for all r > 0. 

By the proof of claim (i) of Theorem [ij for any x G X(V) we can deter- 



mine the coefficients jf(r, to) = 3i{x,r, to) given in (41). Then 

Jj(x,r;t Q ) = jf(x,r,t )F m (x,r) 
are the Jacobi fields along the geodesic 7z(r) that satisfy 
Jj(x,t Q ;t ) = 0, Jj(x, 0;t ) = F m (x,0). 



23 



Let us now consider the set W C M that is a neighborhood of 7x 0ir;o ((0, £o)\ 
C(xq, t}q, to)) and the map : W — >• W given in (11), a point (xo,r ) £ W 



and its small neighborhood V C W so that Xy\v : V — >■ V = Xp-(V) is a dif- 
feomorphism. The inverse of this map defines local coordinates x h-> (sc, r) = 
(X^>|v) _1 (a;) in the set V. We see that the Jacobi fields Jj(x,r) are in fact 
the coordinate vectors for the (x, r) coordinates. Therefore the metric g with 
respect to these coordinates can be recovered by 

#(a^,<3gfc)l(2?,r) = 9mi(x)jf(x,r;to)j l k (x,r;t ). 

Note that here we are in fact varying x, and performing the entire recovery 
of r and s along each geodesic in order to calculate the Jacobi fields along 
that geodesic and can then compute the metric tensor g in the set W in 
the local (x, r) coordinates. Moreover, as we know the coefficients of the 
Jacobi fields Jj(x,r) represented in the parallel frame F m (x,r) along 7$(r), 
we can also find the coefficients of the vectors F m (x, r) in the basis given by 
the Jacobi fields Jj(x,r). Thus we change the local (x,r) coordinates to the 
Fermi coordinates and determine the metric tensor g in the Fermi coordinates 
in some neighborhood W^™^ C W of the set 7x ,r»((0, *o) \ C(x , r] , to)). 
As Wll 

,vo,to ^ S a neighborhood of 7 XQjr?0 ((0, to)\C(xo, rjo, to)) we have not yet 
reconstructed g in the whole neighborhood ofj Xom . To do this, let si > be 
so small that x = 7z ,r?o( — s i) an d Vo = ix ,vo(.~ s i) satisfy Xq £ U and repeat 
the above construction by replacing x by x , r] by rj and t by arbitrary 
to > Si and the spherical surface S by the corresponding surface S . Then, 
we can determine the metric tensor in local coordinates Wl er ™~ . By varying 

x ,r] ,t J jo 

s 1 and t an d using the fact that on a given geodesic the conjugate points of 
a given point form a discrete set, we see that the whole geodesic 7a; 0i7?0 (IR + ) 
can be covered by neighborhoods of the form W/ e l™ . This completes the 
proof of claim (ii) of Theorem p3| 



We finish this section by proving Theorem |4j 

Proof. Let exp x : T X M ->■ M and exp^., : T X ,M' ->■ M' be the exponential 
maps of (M,g) and (M',g'), correspondingly. 

Let p e U and p' G U' be such that p' = $(p) and let 

^ = d^ p ■. T p M Tp/M' 
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Figure 2: Reconstruction procedure in the case of conjugate points. 

be the differential of $ at p. For v = tv° G T P M, \\v°\\ g = 1, t > 0, 
let t„ : T P M — >■ T g M denote the parallel transport along the geodesic 
7py([0, i]), where q = j p , v o(t) in (M,g) and let r^/ : T p iM' — > T q /M' de- 
note the corresponding operation on (M',g'). For v, w G T p M let the curve 
: [0, 2] — )■ M be the broken geodesic, that is defined by 

Hv,w(s) = exp p (sv), for < s < 1, 
M«,t«(s) = exp 9 ((s - l)r„w), for 1 < s < 2, 

where q = exp p (t>). When r = p V)W (2) is the end point of the broken geodesic, 
we denote by t v>w : T p M — >■ T r M the parallel transport of vectors along the 
curve p v>w ([0,2]). For all v G T p M let p(p,v,r) be the function in Theorem 
[3] for the geodesic 7 Pi ^(r), that is, we can determine the Riemannian metric 
in the Fermi coordinates in the tubular neighborhood V XOtVo ^ r that contains 
the ball B^(exp p (v), p(p,v,r)). Let po(v) > be such that po(v) < p(p,v,r) 
and the ball B-^(exp p (v), p (v)) is geodesically convex. Let p' (v') > be 

the corresponding function for (M',g ! ) and p' . Finally, we define f(v) = 
mm(p (v),p' {i{v))). _ _ 

Let v G T p M. When q = exp p (v) and q' = exp' p ,(£(v)), let E v : T p M M 

be the map = ex P(? (r^) and : T p M -> M be = exp'^ «,£')• 

We see that 

o £ o i?- 1 : ^(g, /(„)) -> B^(q f , /(«)) 

is an isometry. In particular, if G T p M are such that ||u>|| 9 < /(f), and 
t = p VtW (2) G M and r' = //^ a m w \(2) G M' are the end points of the broken 
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geodesies, the above implies that the linear isometry 



= T, 



i(v),i(w) 



o£o r~l : T r M -> T r ,M' 



preserves the sectional curvature, i.e., Sec r (£,?7) = Sec r /(£u jM (£), £ V)W (r))). 
Thus from the proof of the Ambrose theorem given in [22] (for the original 
reference, see [2]) it follows that (M,g) and (M',g') have isometric covering 
spaces. 



We begin with the reconstruction of the Riemannian metric g\u if Su 
is given. If x G U then in a sufficiently small neighborhood U' C U of x 
all points z can be connected to x with a geodesic of a given length (travel 
time) contained in U. As a consequence, the distance between x and z can 
be found as 



Indeed, we observe that the infimum is obtained when T,j are the boundaries 
of sufficiently small balls (which are always smooth), Bj^(yj, tj), where yj are 
points on the shortest geodesic connecting x to z. Thus we can determine the 
distance function (y, y') i-> d(y, y') between two arbitrary points in y, y' G U' . 

Now, if r > is small enough and zi,...,z n G U' are disjoint points 
so that d(x,zj) = r, then the function y h-» {d{y,Zj)) T - =1 G W 1 defines local 
coordinates near the point x G U' . So^jn U', we can find the differentiable 
structure inherited from the manifold M. Using the distances beween points 
y,y' G U', we can determine the Riemannian metric in these coordinates in 
U' . But then we can find the Riemannian metric g\u if Su is given. 

For (t, S) G Su and Xo G S, let A^(xo, S, t) be the set consisting of the 
two unit normal vectors of £ at xq. Let A^i(xo,S,t) be the set of those 
Uq G A^(x , S, t) for which the point x has a neighborhood U' C U such that 
SnU' has the representation 



Appendix A: Proof of Proposition [2] 



N 




(43) 



£ntf' = { 7l , )t ,(*); r^G^}, 



(44) 



26 



where y — 7a; ,^o (-t) and Q C is a neighborhood of r/ = -% ,u (-t). 

Note that it is possible for Ni(xo,Y,,t) to contain both normal vectors in 
N(xq, E, t). An example of a case in which this occurs is when M is S 2 and 
E is a subset of the Equator. 

Lemma 7. // £7 and Su are given, we can determine Ni(xo,Yl,t) for any 
(t, E) G 5(7 and so G S. 



Proof. For given (£, E) G «S[/ and x Q G U let Co G iV(x , S, i) be one of the 
two unit normal vectors to E at xo, and let ((x) be a smooth normal vector 
field on E fl U' such that C( s o) = Co- We introduce the notation 



Eg, (a) 



2,±CO) 



»; x e En?/'}; 



E^,(s) will be smooth for s G (— e, e), e > sufficiently small. 



Assume next that Co G -^1(^0 > E, t). Then representation (44) is valid with 
1/ = 7*o,C(«o)( _ *)j and for s ) = 7<n,C(a!)( s ) and ^fo s ) = 7x,cw(s), we have 
7p0m),»;O,*)( _ * - s) = lx£(x)(s ~t - s) = y. Hence, there is a neighborhood 
U" G U of Xq such that for all e > small enough 

(t + s, E+„(s)) G <St/ for all s G (-£, +e). (45) 

Let us consider the following condition 

(C) There exists y' such that I y x ,c(x){+t) = y' for all x G E close to Xq. 

If condition (C) is valid, then both Co and —Co are in N(x , E,i). If condition 
(C) is not valid, then N(xq, E,t) contains the vector Co but not —Co- 
in the case when condition (C) is valid, we see that (45) holds as well as 
the analogous identity with the minus sign, that is, we have 

(t + s, E~„(s)) G Su» for all s G (-£, +e). (46) 

Next, consider the case when the condition (C) is not valid. Our aim is show 
that then (46) can not hold. For this end, let us assume that the condition 
(C) is not valid but we have (46). Then we see that for all s G (— e, e) one of 
the sets 

a+ = K,cw(- S + (* + x g s n [/"}, 

A- = {j x ,c(x)(s - (t + s)); x G E n U"}, 
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would consist of a single point. Now, if A + consisted of a single point then 
this point would satisfy the condition required for the point y' in condition 
(C). As we assumed that condition (C) is not valid, we conclude that A + 
cannot be a single point. If A_ consisted of a single point for all s G (— e, e), 
then for all Xi, x 2 G Ent/" we would have Jxi£(xi)(— t— 2s) = •y X2 ^ X2 )(—t—2s) 
for all s G (— e, e) and hence for all sGl. With s = —t/2 we would see that 



x\ = xi for all x\,x-i G ED U" but that is not possible. Hence equation (46) 
can not be true when the condition (C) is not valid 

Summarizing the above, we can find the set A^i(xq,E, t) using the fact 
that it contains the vector ±( Q if and only if there are U" and e > such 
that (t + s, E^„(s)) G Sjj holds for all s G (— e, e). 

Lemma [7] and the considerations above it prove Proposition [2j 
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